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Abstract 

Origami is the archetype of a structural material with unusual mechanical properties that arise 
almost exclusively from the geometry of its constituent folds and forms the basis for mechanical 
metamaterials with an extreme deformation response. Here we consider a simple periodically folded 
structure Miura-ori, which is composed of identical unit cells of mountain and valley folds with 
four-coordinated ridges, defined completely by 2 angles and 2 lengths. We use the geometrical 
properties of a Miura-ori plate to characterize its elastic response to planar and non-planar piece- 
wise isometric deformations and calculate the two-dimensional stretching and bending response 
of a Miura-ori sheet, and show that the in-plane and out-of-plane Poisson's ratios are equal in 
magnitude, but opposite in sign. Our geometric approach also allows us to solve the inverse 
design problem of determining the geometric parameters that achieve the optimal geometric and 
mechanical response of such structures. 



Folded and pleated structures arise in a variety of natural systems including insect wings 
PP, leaves [2], flower petals [3J, and have also been creatively used by origami artists for aeons 
[1] . More recently, the presence of re-entrant creases in these systems that allows the entire 
structure to fold and unfold simultaneously have also been used in deployable structures 
such as solar sails and foldable maps [5H7j. Complementing these studies, there has been a 
surge of interest in the mathematical properties of these folded structures [U IH [9] , and some 
recent qualitative studies on the physical aspects of origami [TUHT2] . In addition, the ability 
to create them de-novo without a folding template, as a self-organized buckling pattern 
when a stiff skin resting on a soft foundation is subject to biaxial compression [TBHTo] has 
opened up a range of questions associated with their assembly in space and time, and their 
properties as unusual materials. 

Here, we quantify the properties of origami-based 3-dimensional periodically pleated or 
folded structures, focusing on what is perhaps the simplest of these periodically pleated 
structure, the Miura-ori pattern (FigjT^i) which is defined completely in terms of 2 angles 
and 2 lengths. The geometry of its unit cell embodies the basic element in all nontrivial 
pleated structures - the mountain or valley fold, wherein four edges (folds) come together 
at a single vertex, as shown in Fig|Tji. It is parameterized by two dihedral angles 9 G [0, it], 
(3 G [0, it], and one oblique angle a, in a cell of length I, width w, and height h. We treat the 
structure as being made of identical periodic rigid skew plaquettes joined by elastic hinges 
at the ridges. The structure can deploy uniformly in the plane (Figjlp) by having each 
constituent skew plaquette in a unit cell rotate rigidly about the connecting elastic ridges. 
Then the ridge lengths h, I2 and a G [0,7r/2] are constant through folding/unfolding, so 
that we may choose 9 (or equivalently 0) to be the only degree of freedom that completely 
characterizes a Miura-ori cell. The geometry of the unit cell implies that 

/3 = 2sin- 1 (Csin(^/2)), I = 2h(, 

w = 2Z 2 £ and h = l^tsma cos(9/2), 

where the dimensionless width and height are 

f = sinasin(fl/2) and C = cos a(l - £ 2 )~ 1/2 . (2) 

We see that (3, I, w, and h change monotonically as 9 G [0, tc], with (3 G [0, it], I G 2Zi[cos a, 1], 
w G 2/ 2 [0,sina], and h G iifsina, 0]. As a G [0,7r/2], we see that (3 G [9,0], I G [2/i,0], 
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FIG. 1: Geometry of Miura-ori pattern, (a) A Miura-ori plate folded from a letter size paper 
contains 13 by 13 unit cells (along x and y direction respectively), with a = 45° and h = h = Z e - 
The plate dimension is 2L by 2W. (b) In-plane stretching behavior of a Miura-ori plate when pulled 
along the x direction shows its expand in all directions, i.e. it has a negative Poisson's ratio, (c) 
Out-of-plane bending behavior of a Miura-ori plate when a symmetric bending moment is applied 
on boundaries x = ±L shows a saddle shape, consistent with that in this mode of deformation its 
Poisson's ratio is positive, (d) Unit cell of Miura-ori is characterized by two angles a and 6 given 
l\ and li and is symmetric about the central plane passing through O1O2O3. 

w G [0, 21 2 sin(#/2)] and h G [0, Zi]. The geometry of the unit cell implies a number of 
interesting properties associated with the expansion kinematics of a folded Miura-ori sheet, 
particularly in the limit of an orthogonally folds when a = tt/2 (Appendix; A-l), the 
singular case corresponding to the common map fold where the folds are all independent. 
More generally, it is possible to optimize the volume of the folded structure as a function of 
the design variables (A-l). 

From now on, we assume each plaquette is a rhombus, i.e. Zi = Z 2 = Z e , to keep the size 
of the algebraic expressions manageable, although it is a relatively straightforward matter 
to account for variations from this limit. We characterize the planar response of Miura-ori 
in terms of 2 quantities - the Poisson's ratio which is a geometric relation that couples 
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deformations in orthogonal directions, and the stretching rigidity which characterizes its 
planar mechanical stiffness. 

The planar Poisson's ratio is defined as 

_ dw/w 2 



--- dt/t < 3 > 

The reciprocal Poisson's ratio is v lw = Because £ < 1, the in-plane Poisson's ratio 

v wl < (Figj2^i), i.e. Miura-ori is an auxetic material. To obtain the limits on v wU we 
consider the extreme values of a, 9, since u wl monotonically increases in both variables. 
Expansion of ^ shows that v wl \ a -> ~ or 2 \ and thus v wl \e G (— oo, — cot 2 (0/2)], while 
v wl \e^o ~ Q 2 and thus v wl \ a £ (— oo, — cot 2 a]. When (a, 9) = (ir/2,ir), u wl = so that the 
two orthogonal planar directions may be folded or unfolded independently when the folds 
themselves are orthogonal, as in traditional map-folding. Indeed, the fact that this is the 
unique state for which non-parallel folds are independent makes it all the more surprising 
that it is still the way in which maps are folded - since it makes unfolding easy, but folding 
frustrating! Similar arguments can be applied to determine the other geometric Poisson's 
ratios related to height changes, v hl and v wh (A-2.1). 

To calculate the in-plane stiffness of the unit cell, we note that the potential energy of a 
unit cell deformed by a uniaxial force f x in the x direction reads, H = U — Jg f x (dl/d6')d6' , 
where the elastic energy of a unit cell is stored only in the elastic hinges which allow the 
plaquettes to rotate, with U = kl e (9 —6 ) 2 + kl e ((3 — (3 ) 2 , k being the hinge spring constant, 9 
and (3 (= (3 (a, 9 )) being the natural dihedral angles in the undeformed state. The external 
force f x at equilibrium state is obtained by solving the equation 5H/56 = (A-2.2), while 
the stretching rigidity associated with the x direction is given by 

dfx 



4t[(l - Q) 2 + cos 2 a] 



(1 — £q) 2 cos a sin 2 a sin 9 
where £o = £(«,#o) and £ is defined in To understand the limits of K x , we expand 



(A. 8) in the vicinity of the extreme values of a and #o which gives us K x ~ a 2 as a — )■ 0, 
K x ~ (tt/2 - a)- 1 as a -> vr/2, K x ~ 6' 1 as 9 ->■ 0, K x ~ (vr - 9)- 1 as 9 -> vr. We see that 
i^a, has a singularity at (a, 9) = (tt/2,tt). 

We note that K x is not monotonic in either a or 6*o, so that there is an optimal pair of 
these variables for which the stiffness is an extremum. Setting de Q K x \ a = and d a K x \g = 
allows us to determine the optimal design curves, 9 0m (a) (green dotted curve in Fig|2 




(a) (b) 

FIG. 2: In-plane stretching response of a unit cell, (a) Contour plot of Poisson's ratio v wl . v wl 
shows that it monotonically increases with both a and 9. f wl \ a £ [— oo, — cot 2 a], and v wl \o G 
[— oo, — cot 2 (#/2)]. (b) Contour plot of the dimensionless stretching rigidity K x /k. The green 
dotted curve indicates the optimal design angle pairs that correspond to the minima of K x \ a . The 
red dashed curve indicates the optimal design angle pairs that correspond to the minima of K x \g . 
See the text for details. 

and a m (9o) (red dashed curve in Fig|2|3) that correspond to the minimum value of the 
stiffness K x as a function of the underlying geometric parameters denning the unit cell. 
These curves are monotonic, and furthermore 9o m (a) is perpendicular to a = 0, because 
when a — > it is asymptotically approximated by 4(#o m — 7r/2) = a 2 (A-2.3). Similarly, 
a ra (S ) is perpendicular to 9 = 0, because when ^ — it is asymptotically approximated 
by c(a m — a*) = 9%, where c = 4^5 + 5y/E and a* = cos -1 a/ \/E — 2 m 60.9°. Analogous 
arguments allow us to determine the other stretching rigidity K y , which is coupled to K x 
through design angles a and 9 (A-2.2, 2.3). 

To understand the bending response of Miura-ori, we must consider the conditions when 
it is possible to bend a unit cell isometrically, i.e. with only rotations of the plaquettes about 
the hinges. Geometric criteria show that planar folding is the only possible motion using 
rigid rhombus plaquettes in our Miura-ori plates (A-3.1). To enable the bending mode, 
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the minimum model for isometric deformations requires the introduction of 1 additional 
diagonal fold into each plaquette (Fig|3^), either the short fold (e.g. 2 7 ) or the long one 
(e.g. Oi0 8 ). Here, we adopt the short fold as a result of which 4 additional DOFs arise and 
allow both symmetric bending and asymmetric twisting, depending on whether the rotations 
are symmetric or not. 

We see that the out-of-plane bending (FigjTJ^) has Poisson's ratio i/& = —k v /k x > 
[23], where k x and K y are curvatures in the x and y directions. To calculate Uj, in linear 
regime, where the rotations are infinitesimal, we need to first derive the expressions for both 
curvatures. If k x is the curvature in the x direction, it may be expressed as the dihedral angle 
between plane OeO^Og and OaO\0-j (Figj3^i) projected onto the x direction over the unit 
cell length. Similarly, the other curvature component K y may be expressed as the dihedral 
angle between plane O4O5OQ and OjOgOg projected onto the y direction over the unit cell 
width. These are given by 

cos(a/2) S in(e/2) ^ 

= 2i e ^e (02 + 04) ' 



Al e sin(a/2)f 



where 02, 04 are rotation angles about internal folds O7O2 and OsO^ respectively, which are 
positive according to the right-hand rule (A-3.2). We note that although there are a total 
of 5 deformation angles (Figj3^i), both k x and n y depend only on 2 and 4 . This is because 
of the symmetry of deformations about xoz plane; 03 and 05 are functions of 0i and 02 (Eq. 



A. 26 in A), and the case that 0i changes while keeping 2 and 04 being corresponds to the 
planar stretch of a unit cell, so 0i does not contribute to both curvatures. This is consistent 
with our intuition that bending a unit cell requires the bending of plaquettes. The Poisson's 
ratio for bending is thus given by 

^ = -? = -i + r 2 = -^, (6) 

K x 

where the last equality follows from Eqs. ^ and (|5]). If the original plaquettes are allowed 
to fold along the long diagonals instead (e.g. 0$Oi in Fig|3^), the new curvature components 
k x and K y are still given by ^ with a being replaced by tt — a (A-3.3), and 2 , 04 now being 
rotations about axis 8 Oi and O g 2 respectively. Therefore v h = —n y /K x = —v wV This 
result, that the in-plane Poisson's ratio is negative while the out-of-plane Poisson's ratio is 



positive, but has the same magnitude is independent of the mechanical properties of the 
sheet and is a consequence of geometry alone. Although our analysis is limited to the case 
when the deformation involves only small changes in the angles about their natural values, 
it is not as restrictive as it seems, since small changes to the unit cell can still lead to large 
global deformations of the entire sheet. 

Given the bending behavior of a unit cell, we now turn to a complementary perspective 
to derive an effective continuum theory for a Miura-ori plate that consists of many unit cells. 
Our calculations for the unit cell embodied in ^ show that n x j n y is only a function of the 
design angles a and 9, and independent of deformation angles, i.e. one cannot independently 
control k x and K y . Physically, this means that cylindrical deformations are never feasible, 
and locally the unit cell is always bent into a saddle. Mathematically, this means that the 
stiffness matrix of the two-dimensional orthogonal plate [18] is singular, and has rank 1. 
In the continuum limit, this implies a remarkable result: the Miura plate can be described 
completely by a 1-dimensional beam theory instead of a 2-dimensional plate theory. 

To calculate the bending response of a unit cell, we consider the bending stiffness per unit 
width of a single cell in the x direction B x . Although the bending energy is physically stored 
in the 8 discrete folds, it may also be effectively considered as stored in the entire unit cell 
that is effectively bent into a sheet with curvature k x . Equating the two expressions allows 
us to derive B x (A-3.4). In general, B x depends on multiple deformation angles as they 
are not necessarily coupled, although here, we only study the "pure bending" case (A-3.5), 
where a row of unit cells aligned in the x direction undergo the same deformation and the 
stretching is constrained, i.e. <f>i — for all cells, and then <p2 = <Pa must be satisfied. In 
this well-defined case of bending, B x is solely dependent on the design angles, so that 

2" 



BJa,6) =kL 



cot 



16 sin 3 - + I 1 - 



a ( 2 cos a 



1 _ £2)3/2 



(7) 



2/ 2£ 2 cos a sin a cos(#/2) ' 

as shown in Figj3}^, and we have assumed that all the elastic hinges in a cell have the same 
stiffness. 

Just as there are optimum design parameters that allow us to extremize the in-plane 
rigidities, we can also find the optimal design angle pairs that result in the minima of B x , by 
setting dgB x \ a = and d a B x \g = 0. This gives us two curves 8 m (a) and a m {6) respectively 
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(c) (d) 

FIG. 3: Out of plane bending response of a unit cell, (a) The plaquettes deformations about each 
fold are symmetric about the plane O1O2O3, so that the angles 20i, $2, 03, ^4 and 205 correspond 
to rotations about the axes O1O2, 076)25 C^Os, 0$Oz and O3O2 respectively, (b) Numerical 
simulation of the bending of a Miura-ori plate with a = 45° and 9 = 90°. Force dipoles are shown 
by yellow arrows. Color of the folds indicates the value of deformation angles, (c) Contour plot of 
dimensionless bending stiffness B X J (kl e ) corresponding to pure bending of a unit cell. The green 
dotted curve and red dashed curve indicate the optimal design angle pairs that correspond to the 
local minima of B x \ a and B x \g respectively, (d) Contour plot of bending Poisson's ratio. The 
gray scale plot is from the analytic expression [6] and the red curves are extracted from simulation 
results. In our simulations, we use a plate made of 21 by 21 unit cells and vary a from 20° to 70°, 
9 from 30° to 150° both every 10°. 
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shown in Fig. [3j The green dotted curve 8 m (a) starts from (a, 9) ~ (63.0°, 180°), and ends 
at (a, 6) = (90°, 180°). It is asymptotically approximated by 2.2851(a - 1.0995) « (tt- 6 m f 
when a — > 63.0°. The red curve 9 m (a) starts from (a, 0) m (52.3°, 0°), and ends at (a, 9) = 
(90°, 180°), and is asymptotically approximated by 17.7517(a m - 0.9137) « 9 2 when -> 0°. 

The bending stiffness per unit width of a single cell in the y direction B y (A-3.4) is related 
to B x via the expression for bending Poisson's ratio v\ = B x /B y , where z/& is defined in (|6]). 
This immediately implies that optimizing B y is tantamount to extremizing B x . 

The deformation response of a complete Miura-ori plate requires a numerical approach 
because it is impossible to assemble an entire bent plate by periodically aligning unit cells 
with identical bending deformations in both the x and y direction (A-4.1). Our model takes 
the form of a simple triangle-element based discretization of the sheet, in which each edge 
is treated as a linear spring with stiffness inversely proportional to its rest length. Each 
pair of adjacent triangles is assigned an elastic hinge with a bending energy quadratic in 
its deviation from an initial rest angle that is chosen to reflect the natural shape of the 
Miur-ori plate. We compute the elastic stretching forces and bending torques in a deformed 
mesh [THl ED], assigning a stretching stiffness that is six orders of magnitude larger than the 
bending stiffness of the adjacent facets, so that we may deform the mesh nearly isometrically 
(A-4.2). When our numerical model of a Miura-ori plate is bent by applied force dipoles 
along its left-right boundaries, it deforms into a saddle (Fig|3}D). In this state, asymmetric 
inhomogeneous twisting arises in most unit cells; indeed this is the reason for the failure of 
averaging for this problem since different unit cells deform differently. This is in contrast 
with the in-plane case, where the deformations of the unit cell are afnnely related to those 
of the entire plate. 

To compare the predictions for the bending Poisson's ratio i/f, of the one-dimensional 
beam theory with those determined using our simulations, in Fig|3]i we plot v\> from Q 
(the gray scale contour plot) based on a unit cell and v h extracted at the center of the bent 
Miura-ori plate from simulations (the red curves). We see that these two results agree very 
well, because the unit cell in the center of the plate does have a symmetry plane so that 
only symmetric bending and in-plane stretching modes are activated, consistent with the 
assumptions underlying (|6]). (A-4.2.) 

Our physical analysis of the properties of these folded structures, mechanical metama- 
terials that might be named Orikozo, from the Japanese for Folded Matter are rooted in 
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geometry of the unit cell as characterized by a pair of design angles a and 6 together with 
its symmetry and the constraint of isometric deformations. It leads to simple expressions for 
the linearized planar stretching rigidities K x , K y , and non-planar bending rigidities B x and 
B y . Furthermore, we find that the in-plane Poisson's ratio v wl < 0, while the out-of-plane 
bending Poisson ration vj, > 0, an unusual combination that is not seen in simple materials, 
satisfying the general relation i.e. v wl = —Uj,; a consequence of geometry alone. Our analysis 
also allows us to pose and solve a series of design problems to find the optimal designs of 
the unit cell that lead to extrema of stretching and bending rigidities as well as contrac- 
tion/expansion ratios of the system. This paves the way for the use of optimally designed 
Miura-ori patterns in such passive settings as three-dimensional nanostructure fabrication 
[2~T] . and raises the possibility of optimal control of actuated origami-based materials in soft 
robotics [22] and elsewhere using the simple geometrical mechanics approaches that we have 
introduced here. 

We thank the Wood lab for help with laser cutting to build the paper Miura-ori plates 
shown in Figure [TJ and the Wyss Institute and the Kavli Institute for support, and Tadashi 
Tokieda for many discussions and the suggestion that these materials be dubbed Orikozo. 
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APPENDIX 



1. GEOMETRY AND KINEMATICS 

Before we discuss the coupled deformations of the plate embodied functionally as {3(a, 9), 
we investigate the case when a = tx/2 corresponding to an orthogonally folded map that 
can only be completely unfolded first in one direction and then another, without bending 
or stretching the sheet except along the hinges. Indeed, when a = tx/2 and 9 ^ tt, Eq. 
reduces to /3 = 0, I = and h = l\, the singular limit when Miura-ori patterned sheets 
can not be unfolded with a single diagonal pull. Close to this limiting case, when the folds 
are almost orthogonal, the Miura-ori pattern can remain almost completely folded in the x 
direction (/3 changes only by a small amount) while unfolds in the y direction as 9 is varied 
over a large range, only to expand suddenly in the x direction at the last moment. This 
observation can be explained by expanding Eq. Q asymptotically as a — > tt/2 and 9 — > tt, 
which yields « tt - e/8, lmh(2- (e/5) 2 /4), w « Z 2 (2 - 5 2 - e 2 /4) and h w ke/{28), where 
5 = tt/2 — a and e = ix — 9. Thus, we see that for any fixed small constant 5, only when 
e < 5, do we find that /? — > tt, I — > 2li and h — > 0, leading to a sharp transition in the narrow 



neighborhood (~ 5) of 9 = n as a — )■ 7r/2 (Fig.A.li), consistent with our observations 



More generally, we start by considering the volumetric packing of Miura-ori characterized 
by the effective volume of a unit cell V = I xw x h — 2Z 2 Z 2 C 2 si n & sm a t an a i which vanishes 
when 9 = 0, tt. To determine the conditions when the volume is at an extremum for a fixed 
in-plane angle a, we set deV\ a = and find that the maximum volume 

2 2 _, / cos 2a — 1 \ 

Vmax \a = 2/ 1 / 2 sm a at m = cos — — , (A.l) 

\cos 2a + 3 J 

shown as a red dashed line in Fig jA.lb . Similarly, for a given dihedral angle 9, we may ask 
when the volume is extremized as a function of a? Using the condition d a V\g = shows 
that the maximum volume is given by 

Al\l 2 cos a m (VB + 4 cos 9 - 3) cot 2 (9/2) sin 

Kiai = /r , = ~ (A. 2) 

V5 + 4 cos 6» - 3 - 2 cos 6» 

at 

-i 



a m = cos 



z + cos 9 - V5 + 4 cos e) /{cos 9 - 1) 
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FIG. A.l: Geometry of the unit cell as a function of a and 8. (a) The folding angle /3 increases as 
8 increases and decreases as a increases. The transition becomes sharper as a ~ vr/2, and when 
a = vr/2, (3 = independent of 8, i.e. the unfolding (folding) of folded (unfolded) of maps with iV 
orthogonal folds has 2 N decoupled possibilities, (b) Effective dimensionless volume V/{2l\l2). The 
green dotted curve 8 m {a) indicates the optimal design angle pairs that correspond to the maximum 
V\ a . The red dashed curve a m {8) indicates the optimal design angle pairs that correspond to the 
maximum V\q. (c) Contour plot of the dimensionless stretching rigidity K y /k. K y \ a is monotonic 
in 8q. The green dotted curve indicates the design angle pairs that correspond to the minima of 
K y \o . The red dashed curve indicates the design angle pairs that correspond to the maxima of 
K x \q . See the text for details. 



shown as a red dashed line in FigjATp. These relations for the maximum volume as a 
function of the two angles that characterize the Miura-ori allow us to manipulate the con- 
figurations for the lowest density in such applications as packaging for the best protection. 
In the following sections, we assume each plaquette is a rhombus, i.e. h = h = l e , to keep 
the size of the algebraic expressions manageable, although it is a relatively straightforward 
matter to account for variations from this limit. 
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2. IN-PLANE STRETCHING RESPONSE OF A MIURA-ORI PLATE 



2.1 Poisson's ratio related to height changes 

Poisson's ratios related to height changes, u hl and u wl read 

i dh/h o 2 ^ 
v,,—v = — — = cot a sec -, 

2 fA3) 

i/. = z/^ 1 = — — V- = C tan 2 -. 
which are both positive, and monotonically increasing with 9 and a. Expansion of v in Eq. 



(A. 3) shows that u hl \g^. n ~ (vr — #)~ 2 and thus v M \ a G [cot 2 a, oo), while v hl \ a ^o ~ and 
thus u hl \g G (oo, 0]. Similarly, expansion of u hw in Eq. (A. 3) shows that v hw \g^ ~ (n — 6)~ 2 
and thus v hw \ a G [0,oo), while G [tan 2 (0/2), 0]. Finally, it is worth pointing out that 

has a singularity at (a, 6) = (tt/2,tt). 

2.2 Stretching stiffness K x and 

Here we derive the expressions for stretching stiffness K x and K y . 

The expression for the potential energy of a unit cell deformed by a uniaxial force f x in 
the x direction is given by 

r e dl 

H = U - I f*dV d6 '> ( A - 4 ) 
where the unit cell length I is defined in Eq. The elastic energy of a unit cell U is stored 
only in the elastic hinges which allow the plaquettes to rotate, and is given by 

U = kl e (9-6 ) 2 + kl e ((3-p ) 2 , (A.5) 

where k is the hinge spring constant, and 9 and f3 (= {3(a, 9 )) are the natural dihedral 
angles in the undeformed state. The external force f x at equilibrium state is obtained using 
the condition that the first variation 5H/58 = 0, which reads 

_ dU/dO _ nu (0-6 o ) + (p-i3 o )zj(a,6) 
fx ~ dl/dO ~ Zk v (a,9) ' [ } 



where U is defined in Eq. ( |A.5 ), I is defined in Eq. ([T]), and in addition 



cos a , . cos a sin 2 a sin 9 

^(a,9) = ^—^ and r)(a, 9) = _ 2 . (A.7) 
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The stretching rigidity associated with the x direction is thus given by 



K x (a,9 ) = 



d9 



4k- 



1 - + cos 2 a 



(A.8) 



[1 — £q) 2 cos a sin 2 a sin #0 

where £ = £(a,9 Q ). 

Similarly, the uniaxial force in the y direction in a unit cell at equilibrium is 

= dU/dO = 2k {9-0 Q ) + {P-p )w{a,0) 
y dw/d8 sinacos(#/2) 

where w is defined in Eq. Q and w is defined in Eq. (A.7). The stretching rigidity in y 
direction is thus given by 



(A.9) 



K y (a,9 ) ee 



dfy 
dQ 



2k- 



(l-e f + cos 2 a 
1 -^) 2 sinacos(#o/2y 



(A.10) 



of which the contour plot is show in Fig. A.l 



2.3 Asymptotic cases for optimal design angles 

The expressions in Section 2.2 allow us to derive in detail all the asymptotic cases associ- 
ated with the optimal pairs of design angles which correspond to the extrema of stretching 
rigidities K x and K y . For simplicity, we use (a, 6) instead of (a, 60) to represent the design 
angle pairs when the unit cell is at rest. 

1. Expanding dgK x in the neighborhood of a = yields 

o T ^ 1 8 cot 6 esc 9 2 , m 2 „x _ 1 . 2 . ,a i,\ 

MkU-x] = o ((3 + cos9)csc 2 9) +0(a 2 ). (A.ll) 



As a — 0, 6* — 7- 7r/2 to prevent a divergence. Continuing to expand Eq. (A.ll) in the 
neighborhood of 9 = 71 /2 and keeping the first two terms yields 

d e K x \ e ^ /2 = 4(6 - vr/2) = a 2 . (A.12) 

Therefore in the contour plot of K x (Fig.3b in the main text), the greed dotted curve is 
approximated by 4(9 — 7r/2) = a 2 in the neighborhood of a = 0, and is perpendicular 
to a = as 9 is quadratic in a. 
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2. Expanding d a K x in the neighborhood of 9 = yields 

_ . [11 + 20 cos(2a) + cos(4a)l esc 3 a sec 2 a 1 ,, . 
d a K x \ ^ o = - [ - oa 77^ {t 290 + 173cos(2a) 

+46 cos(4a) + 3 cos(6a)] esc 3 a sec 2 a} 9 + 0(9 2 ). 

(A.13) 



The numerator of the leading order in Eq. (A.13) has to vanish as 9 — > to keep the 



result finite, which results in a unique solution a* = cos 1 ^y/ y/E — 2^ in the domain 



a G (0, 7r/2). Continuing to expand Eq. (A.13) in the neighborhood of a = a* and 
only keeping the first two terms yields 



d a K x = 0\ a ^ a * => 4^5(1 + VE)(a - a*) = 9 2 . (A.14) 

so the red dashed curve in the contour plot of K x (Fig.3b in the main text) is perpen- 
dicular to 9 = 0. 

3. Similarly, Expansion of d a K y near 9 = n yields 

_ i [— 1 + 16 cos(2a) + cos(4a)l esc 2 a sec 3 a 1 r . 

d a K y \ e ^ = [ - + 77^ 638 - 737cos 2a + 

2(9 — tt) 192 

162cos(4a) + cos(6a)] esc 2 a sec 5 a(9 - tt) + 0[(9 - tt) 3 ]. 

(A.15) 

Allowing for a well behaved limit at leading order as 9 — > n requires —1 + 16 cos(2o;) + 
cos(4a) = and yields a* = cos -1 (yj (y/l7 — 3)/2^ as the unique solution when a 



is an acute. Again expanding Eq. (A.15) in the neighborhood of 9 = tt, and only 
keeping the first two terms yields 



d a K y \ a ^ a , = 2 v/l + Vl7(a m - a*) = (tt - 9f. (A.16) 



So the green dotted curve in the contour plot of K y (Fig. A.l:) is approximated by 
2\/l + y/l7(a m — a*) = (tt — 9) 2 near a = a*, and is perpendicular to 9 = tt. The 
point where the green curve ends satisfies the condition 

d a K y = and d a (d a K y ) = (A. 17) 

and numerical calculation gives us the coordinates of this critical point as 

9 = 2.39509, and a = 1.00626. (A.18) 

The red dashed curve (Fig. A.lfc ) starting at this point shows a collection of optimal 
design angle pairs (at, 9) where K y \g is locally maximal. 
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3. OUT-OF-PLANE BENDING RESPONSE OF A MIURA-ORI PLATE 



3.1 Minimum model for isometric bending 

Here we show that planar folding is the only geometrically possible motion under the 
assumption that the unit cell deforms isometrically, i.e. with only rotations of the rhombus 
plaquettes about the hinges. To enable the out-of-plane bending mode, the minimum model 
for isometric deformations requires the introduction of 1 additional diagonal fold into each 
plaquette, and this follows from the explanation below. 



Suppose the plane Oi0 2 5 4 (Fig, A. 2 1) is fixed to eliminate all rigid motions, for any 
dihedral angle 9, the orientation of plane O1O2OSO7 is determined. However, the other two 
rhombi O2O5OQO3 and O2O3O9O8 are free to rotate about axis O2O5 and O2OS respectively 



and sweep out two cones which intersect at O2O3 and O2O3. Fig A. 2 a, shows the two possible 
configurations of a unit cell determined from the two intersections, the yellow part being the 
red part that has been flipped about a plane of symmetry. The unit cell in red is the only 
nontrivial Miura pattern, so that for any given 9, there is a unique configuration of the unit 
cell corresponding to it. Any continuous change in 9 results in the unit cell being expanded or 
folded but remaining planar, in which case, Oi, O4, O7, O3, 0§ and Og also remain coplanar. 
In order to enable the bending mode of the unit cell, the planarity of each plaquette must 
be violated. In the limit where the plaquette thickness tcl the stretching rigidity (~ t) is 
much larger than the bending rigidity (~ t 3 ), with t being the thickness of a plaquette, while 
the energy required to bend a strip of ridge is 5 times of that required to stretch it according 
to the asymptotic analysis of the Foppl — von Karman equations [23]. Therefore, the rigid 
ridge/fold is an excellent approximation for out-of-plane bending when t <C 1. Then, to get 
a bent shape in a unit cell and thence in a Miura-ori plate, we must introduce an additional 



fold into each rhombus to divide it into two elastically hinged triangles (Fig, A.2p ). As a 



result, 4 additional degrees of freedom are introduced in each unit cell. The deformed state 



can either be symmetrical about the plane O1O2O3 (Fig, A. 2 2) corresponding to a bending 



mode, or unsymmetrical corresponding to a twisting mode. Here, we are only interested in 



the bent state, in which the rotation angle 02 about the axis O2O4, and 04 about the axis 



O3O5, are the same as the rotations about 7 2 and 8 3 respectively. The rotation angles 



about the axis 0iO 2 , 3 2 and 2 5 are 20! , 20 5 and 3 respectively. (— > indicates the 
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(a) (b) (c) 

FIG. A. 2: Bending of a unit cell, (a) The two configurations of a unit cell for any given 9 if 
each plaquette is a rigid rhombus. The only possible motion is in-plane stretching. The yellow 
plaquettes illustrate the trivial configuration of two rigid plaquettes and the red ones show the 
typical configuration of a Miura-ori unit cell, (b) The undeformed state. An additional fold along 
the short diagonal is introduced to divide each rhombus into 2 elastically hinged triangles, (c) 
Symmetrically bent state. The bending angles around axis O2O4 and O3O5 are the same as those 
around O7O2 and 0$Os respectively. 

direction.) 

3.2 Curvatures and the bending Poisson's ratio when short folds are introduced 

Here we derive expressions for the coordinates of every vertex of the unit cell after bending 
in the linear deformation regime, from which curvatures in the two principle directions k x , 
K y and the bending Poisson's ratio i/j, = —K, y /n x can be calculated. 

To do so, we first need to know the transformation matrix associated with rotation about 
an arbitrary axis. The rotation axis is denned by a point {a, b, c} that it goes through 
and a direction vector < u,v,w >, where u, v, w are directional cosines. Suppose a point 
{xq, yo, Zq} rotates about this axis by an infinitesimal small angle u (u < 1), and reaches the 
new position {x,y, z}. Keeping only the leading order terms of the transformation matrix, 



18 



we find that the new position {x, y, z} is given by 



x = Xq + (— cv + bw — wy + vz )ui, 

V — VO + ( cu ~ aw + wx — M2o)w, 

z = z + (—bu + av — vxq + uyo)uj. 



(A.19) 



Given Eq. (A.19), we are ready to calculate the coordinates of all vertices in the bent 



sate. Assuming that the origin is at 2 , in the undeformed unit cell, edge 0\0<i is fixed in 
xoz plane to eliminate rigid motions. Each fold deforms linearly by angle 20i, 2 , 03, 04 
and 205 (see Fig. A.2fc ) around corresponding axes respectively. The coordinates of 0\ and 
O2 are 

sa sin a cos(0/2) 

!-^ 2 V 1 "^ 2 (A.20) 

O2* = 0, o 2y = 0, 22 = 0. 



o 



Ix 



o 3x 

3 y 



The coordinates of O3 after bending are 

cos a cos(a/2) sin a sin0 



+ 



sin a sin 



1 - sin 2 a sin 2 (6/2) ^3 - cos(2a)(cos0 - 1) + cos0 2 y/l - sin 2 a sin 2 (6/2) 
4cos(0/2) sin(2a) 
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3 + cos(2a) + 2 cos 6 sin 2 a 

+ cos(6 l /2) sin(a)0 3 , 

cos(6 l /2) sin a 2 cos a cos(a/2) sin(6 l /2) 



csc(0/2)[-sina; + sin(2a) + sin 3 a sin 2 (6 l /2)] sin0 
1 + [3 + cos(2a) + 2 cos(0) sin 2 a] sin(a/2) ^ 2 



+ 



cosasinasm(0/2) 

- sin 2 a sin 2 (6/2) ' ^3 - cos(2a) (cos - 1) + cos 2 ^1 - sin 2 a sin 2 (0/2) ' 

(A.21) 

The coordinates of O4 after bending are 
cos a + sin 2 a sin 2 (6/2) - 1 



sin 2 a sin 

V^l - sin 2 a sin 2 (0/2) ^3 - cos(2a) (cos 6 - 1) + cos 6 U 

O/iy = sin a sin(0/2) — cos(0/2) sino;0i, 

cos(0/2) sin a 2cosasinQ;sin(0/2) 



(A.22) 



O 



Az 



y/l - sin 2 a sin 2 (0/2) ^3 - cos(2a)(cos0 - 1) + cos0 
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The coordinates of 5 after bending are 

^ r. ~o . 9//1 7T7 sin 2 a sin 

5x = - VI - sm 2 asm 2 (0/2) - - 4 
v a/3 — cos(2a)(cos0 — 1) + cos 6* 

sin 2 a sin 

+ 



2a/3 - cos(2a)(cos0 - 1) + cos 6* sin(a/2) 

^ • • m m\ { r\ i c\\ ■ , cos(0/2) sina , 

C 5 „ = sina sm(6'/2) — cos(6'/2) sin a®! H — — - — <f> 2 , 

2sm(a/2) 

^ 2cosasinasin(0/2) , cos a sin a sin(0/2) 

5z = w 01 + 02- 

a/3 — cos(2a)(cos0 — 1) + cos0 a/3 — cos(2a)(cos 9 — 1) + cos0 sin(a/2) 

(A.23) 

The coordinates of Oq after bending are 

sin 2 a sin 2 (0/2) — cos a — 1 sin 2 a sin 



a 



6j; 



a/1 - sin 2 a sin 2 (9/2) a/3 - cos(2a) (cos 9 - 1) + cos 9 
sin 2 a sin sinasin0cos(a/2) 



2 a/1 - sin 2 a sin 2 (9/2) a/3 - cos(2a)(cos 9 - 1) + cos0 

^ . . ,„ 4cos(0/2)sina[sin 2 asin 2 (0/2) - 1 - 2cosal , 

O ey = sin a sin 0/2 + 1 ' ' 1 1 V- 2 Vi 

3 + cos(2a) + 2 cos 9 sin a 

8 cosa cos(0/2) cos(a/2) , ■ , / , ✓» 

+ ^ , y J a - 2 & + cos(0/2) sin a0 3 - cos(0/2) cos(a/2)0 4 , (A.24) 

6 + cos(2a) + 2 cos 9 sin a 

cos(0/2)sina 2 cos a sina sin (0/2) 

Cfi z = — — , = — — , 01 

\/l — sin 2 a sin 2 (0/2) A/3-cos(2a)(cos0-l) + cos0 

csc(a/2) sin(2a) sin(0/2) cos a sina sin(0/2) 

a/3 - cos(2a)(cos0 - 1) + cos0 2 a/1 - sin 2 a sin 2 (0/2) 3 

2 cos a cos(a/2) sin(0/2) 

H , ^=0 4 . 

a/3 — cos(2a)(cos0 — 1) + cos0 



(A.25) 



The coordinates of O7, 0$ and Og after bending are 

{0 7x , 7y , 7z } = {O iX) — 4y , 4z }, {0 8x , 8y , 8z } = {0 5x , —0 5y , 5z }, 

{C>9 X , Og y , Og z } = {Oq x , —Oq v , Oq z }. 

Due to symmetry, 3 must lie in the xoz plane after bending, so 3y = 0, from which 3 
and 05 can be expressed as a function of 0i and 2 , 

8cosa 1 /a\ / 8cosa \ 

~ 3 + cos(2a) + 2cos0sin 2 a + 2 CSC V 2 J V ~ 3 + cos(2a) + 2 cos sin 2 a / 



1 /a 
>i — - esc 



(A.26) 
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The curvature of the unit cell in the x direction is defined as the dihedral angle formed 
by rotating plane 4 0i0 7 to plane 6 3 9 projected onto the x direction over the unit 
length /. The sign of the angle follows the right-hand rule about the y axis. The dihedral 
angle between plane O4O1O7 and plane xoy is 

4z -0 lz 4 cos a sin a sin (9/2) /*o^ 

SZ417 = — . = ; — k 01, (A. 27 

V /T3^2 3 + cos(2a) + 2cos#sin 2 a v ; 

and the dihedral angle between plane 3 6 9 and plane xoy is 

n = Oez - Ozz = 2[cos(a/2) + cos(3a/2)][0 2 + 4 - 20! sin(a/2)] sin(0/2) 
639 y/l^ 3 + cos(2a) + 2 cos sin 2 a ' [ ' ' 

The curvature k x hence is 

^639-^417 (02 + 4 )cos(a/2)sm(0/2) 
«. = 1 = • (A.29) 

The curvature of the unit cell in the y direction is defined as the dihedral angle between 
plane O4O5O6 and OjO^O^ projected onto the y direction over the unit cell width w, which 
is expressed as 

205^-0^-03^ 1,, f a.\ re 



S = h * = -4(^2 + 4 ) csc ( - ) csc ot csc ( - ) Vl - e- (A.30) 



From Eq. (A.29) and Eq. (A.30), we can calculate the bending Poisson ratio, which is 
simplified to 

u b = -^L = -1 + csc 2 a csc 2 ( - ) . (A.31) 
k x \2 J 

3.3 Curvatures and the bending Poisson's ratio when long folds are introduced 



In Fig]A.2[ if we introduce the additional fold along the long diagonal, e.g. O1O5, instead 



of the short one, the unit cell can be bent too. In this case, 02 and 4 are bending angles 



around axis O1O5 and O2O6 respectively. Oi, O2 do not change as they are fixed, and 
coordinates of O3 after bending are 

„ cos a sin 2 a sin , sinasin(a;/2) sin0 

O,,,- = == + 



1 - sin 2 a sin 2 (0/2) 20 - sin 2 a sin 2 (9/2) ^3 - cos(2a)(cos 9 - 1) + cos0 

4cos(0/2)sin(2a) , f 9\ . . f9\ . /a\ , 

3 + cos(2a) + 2 cos 9 sin a \2y \2y V2/ 

cos(0/2) sin(a) cosa;smasm(0/2) 2cosasin(a/2) sin(0/2) 

03^ = 03 H 0, 

VI -sin 2 a sin 2 (0/2) V 1 

— sin 2 a sin 2 (0/2) y 3 — cos(2a)(cos0 — 1) + cos 6* 

(A.32) 
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The coordinates of 4 after bending are 

cos a — 1 + sin 2 a sin 2 (9/2) sin 2 a sin 



Oa x — 



a/1 - sin 2 a sin 2 (9/2) a/3 - cos(2a)(cos0 - 1) + cos0 

O^y = sin a sin(0/2) — cos(0/2) sino:0i + cos(0/2) sin(a/2)0 2 , 
cos(6 l /2) sin a 2cosasin(0/2) sin(ai/2) 



0i-^sec(|) 



O. 



4z 



a/I -sin 2 a sin 2 (0/2) a/3 - cos(2a) (cos 9 - 1) + cos 9 
The coordinates of 5 after bending are 



2 COS I ^ j 01 — 02 



(A.33) 



= — W 1 — sin 2 a sin 2 



sin asm i 



2/ a/3 - cos(2a) (cos 9 - 1) + cos 9 ' 
5j/ = sin a sin(6 l /2) — cos(0/2) sin «0i, 



(A.34) 



52 = 



sin(2a) sin(0/2) 



= 01- 



a/3 — cos(2a)(cos — 1) + cos 
The coordinates of 6 after bending are 

sin 2 a sin 2 (0/2) — cos a — 1 sin 2 a sin 9 



6x 



01 + 



sin 2 a sin 9 



a/1 - sin 2 a sin 2 (0/2) a/3 - cos(2a)(cos - 1) + cos0 2a/1 - sin 2 a sin 2 (0/2) 



Oq v = sin a sin 



0\ 4cos(0/2)sina [-1 - 2 cos a + sin 2 a sin 2 (0/2)1 /0\ 

t: H ; — ■ ^ o 0i + cos - sino;03, 

2) 3 + cos(2a) + 2cos0sin 2 a Y \ 
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cos(0/2) sin ck 



sin(2a)sin(0/2) 



=0i 



cos a; sin a sin(0/2) 



a/1 - sin 2 a sin 2 (0/2) a/3 - cos(2a) (cos 9 - 1) + cos 9 a/1 - sin 2 a sin 2 (0/2) 

(A.35) 

Using the same idea for the long fold case as we did for the short fold, we can also 
calculate the curvatures in the two principal directions and find that 

^639 - ^417 2[sin(a/2) - sin(3a/2)] sin(0/2) , . sin (a/2) sin (0/2) 

K.x = : = —TZ 7^- : „ . o ~. (02 + 04) = ~, (02 + 04), 



I 



[3 + cos(2a) + 2 cos sin 2 a]l 



2l e ^l^e 



(A.36) 



while 



ZO^y — O^y — O^y 

hw 



sin a sin 



2 (0/2) 



2 cos(a/2)w 



(02 + 04 



4/ e cos (a/2)f 



(02 + 04 



Therefore the bending Poisson ratio is 



V h = -- 



K,, 



— 1 + csc a esc 



(A.37) 



(A.38) 



which is the same as that of the case when the short folds are introduced. 
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3.4 Bending stiffness B x and B y 

We are now in a position to derive expressions for the bending stiffness B x and B y . On 
one hand, the bending energy is physically stored in the 8 discrete folds, which can be 
expressed as l/2kl e [A4>\ + 4sin(a/2)</>2 + 20| + 4sin(a/2)0| + 40|]. On the other hand from 
a continuum point of view, the energy may also be effectively considered as stored in the 
entire unit cell that is bent into the curvature k x , which can be expressed as 1/2B x wIk x . 
Equating the two expressions for the same energy, we can write B x as 



2 



= m + 4sin(f )fi + 24>l + 4sin(f )0l + 40 5 



Similarly, the bending stiffness per unit width of a single cell in the y direction is 



if + 4 sin(f )<fe + 20|_+ 4 sin(f )0l + A ' 1 
win} 



B y (a,9) = kl e -^ V2/ ^ • ( A - 4 °) 



3.5 Pure bending 

Finally, we explain the "pure bending" situation in the main text, borrowing ideas from 
notions of the pure bending of a beam where curvature is constant. If we demand that a row 



of unit cells aligned in the x direction (e.g. the cell C\ and in Fig. A. 3) undergo exactly 
the same deformation, this results in 02 = 4>a- Furthermore, in this limit, the stretching 
mode is constrained, so that <pi — for all cells. For this well defined bending deformation, 
the bending stiffness depends only on the design angles, not on the deformation angles as 



shown in Eq. (A. 39) and Eq. (A.40) 



4. NUMERICAL SIMULATIONS OF THE BENDING RESPONSE OF A MIURA- 
ORI PLATE 

4.1 Homogeneous deformation in bent plate is impossible 

Here we explain why it is impossible to assemble an entire bent plate by periodically 
aligning unit cells with identical bending deformation in both the x and y direction. 



In Fig, A. 3 , the 4 unit cells C\, C2, C3 and C4 have identical bending deformations: G\ 
and C 2 align perfectly in the x direction, which requires that ZO^OiOj = Z0 6 3 0g = 
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ZO11O13O15. G\ and C3, C% and C4 align perfectly in the y direction respectively, which is 
automatically satisfied by the symmetry of the unit cell. Now the question becomes whether 
the unit cell C3 and C4 can align together? The answer is no. The reasoning is as follows. 
O3 and 3 are symmetric about plane OqOuOi^, while O3 and 0' 3 ' are symmetric about 
plane O^O^Oq. However plane O6O12O13 and plane O^O^Oq are not coplanar unless all the 
deformation angles about the internal folds are zero, which is violated by bending. 0' 3 and 



3 thus do not coincide. In fact 3 = 3 if and only if 3y 



O 



6y , which requires 



from Eq. (A.22), Eq. (A.23), Eq. (A.24) and Eq. (A.26). This is the in-plane 



stretching mode instead of the bending mode. In conclusion, in the bent Miura-ori plate, 
the deformation must be inhomogeneous. 




FIG. A. 3: 4 unit cells with identical bending deformation cannot be aligned together to form an 
entire plate. See the text for details. 



4.2 Simulation model 

In this subsection, we explain the bending model and the strategies used to bend the 
Miura-ori plate. 

We endow these triangulated meshes with elastic stretching and bending modes to cap- 
ture the in-plane and out-of-plane deformation of thin sheets. The stretching mode simply 
treats each edge in the mesh as a linear spring, all edges having the same stretching stiff- 
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(a) (b) 

FIG. A. 4: Simulation model, (a) A single bending adjacency. The vectors Uj illustrate the purely 
geometric bending mode and Ni and N 2 are the weighted normals of the adjacent triangles, (b) 
The left-right bending strategy is shown in yellow and the up-down bending strategy is shown in 
green. Each arrow represents a force applied to its incident vertex. Left-right force directions bisect 
the yellow adjacencies and are perpendicular to the shared edge and up-down force directions are 
normal to the plane spanned by each pair of green edges. 

ness. Accordingly, the magnitude of the restorative elastic forces applied to each node in a 
deformed edge with rest length x and stretching stiffness k is given by ^(x' — Xq) and the 
energy contained in a deformed edge is given by 

^(x'-xo) 2 . (A.41) 

The Xq term in denominator of the stretching mode ensures mesh-independence. The bend- 
ing mode is characterized in terms of four vectors U\, u 2 , u 3 an d U4, each of which is 
applied to a node in a pair of adjacent triangles. Defining the weighted normal vectors 
Ni = (xi — £3) x (xi — X4) and N 2 = (x 2 — 24) x (x 2 — X3) and the shared edge E = 24 — X3, 
we may write 



"1 = \E\j^ (A.42) 



(xt-x^-E A r i (x 2 -x 4 )-E N 2 (\aa\ 



LEI \NA 2 \E\ \N. 



2 



(xi - x 3 ) • £ A4 (x 2 — x 2 ) ■ E N : 



u4 = __^_^___^^ (A45) 
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The relative magnitudes of these vectors constitute a pure geometric bending mode for a 
pair of adjacent triangles. For pairs of adjacent triangles that do not straddle the fold line, 
the force on each vertex is given by 

Fi = h( 6 - - |K, (A.46) 

where k b is the bending stiffness and 9 is the angle between Ni and N 2 that makes each Ui 
a restorative force. For pairs of adjacent triangles that straddle folds, 9 is non-zero and 
shifts the rest angle of the adjacency to a non-planar configuration. The bending energy 
contained in a pair of adjacent triangles is given by 

E b = k b f d -- %9, (A.47) 



with a precise form of 



E b = hC- - |) 2 , (A.48) 



which is quadratic in 9 for 9 ~ 9 . 

We introduce viscous damping so that the simulation eventually comes to rest. Damping 
forces are computed at every vertex with different coefficients for each oscillatory mode, 
bending and stretching. We distinguish between these two modes by projecting the velocities 
of the vertices in an adjacency onto the bending mode, and the velocities of the vertices in 
an edge onto the stretching mode. 

We use the Velocity Verlet numerical integration method to update the positions and 
velocities of the vertices based on the forces from the bending and stretching model and the 
external forces from our bending strategies. At any time t + At during the simulation we 
can approximate the position x(t + At) and the velocity x(t + At) of a vertex as 

x(t + At) = x(t) + x(t) At + ^x(t) At 2 , 

«t + A^m + m+ f +&t, At. (A ' 49) 

A single position, velocity and accleration update follows a simple algorithm. 

• Compute x(t + At) 

• Compute x(t + At) using x(t + At) for stretching and bending forces and x(t) for 
damping forces 
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• Compute x(t + At) 

Note that this algorithm staggers the effects of damping on the simulation by At. 

In simulation, the Miura-ori plate is made of 21 by 21 unit cells, 21 being the number 
of unit cells in one direction, a varies from 20° to 70°, and 6 varies from 30° to 150°, 
both every 10°. We design two bending strategies, each of which corresponds to a pair of 
opposite boundaries. The left-right bending strategy identifies the adjacencies with 2 3 
shared edges on left boundary unit cells and O1O2 shared edges on right boundary unit 
cells (highlighted in yellow in Fig, A.4b ). For each of these adjcencies we apply equal and 



opposite forces to the vertices on their shared edge, the directions of which are determined 
to lie in the bisecting plane of O1O2O4 and O1O2O7 (left boundary unit cells) and O2O3O5 
and 2 3 8 (right boundary unit cells) and perpendicular to the shared edge. The up-down 
bending strategy identifies the top edges of each unit cell on the up and down boundaries 
of the pattern (shown in green in Fig A. 4b ). Each unit cell has one such pair of edges and 



we apply equal and opposite forces to the not-shared vertices in this pair, the directions 
of which are normal to the plane spanned by the pair of edges. We take out the 11th row 
and 11th column of vertices on the top surface as two sets of points to locally interpolate 
the curvature near the center of the plate in x and y direction respectively. The largest 
difference of v h for the same design angle pairs a and 6 between both B.Cs applied is less 
than 0.5%. 

By applying the bending strategies described above, we are able to generate deformed 



Miura-ori plates in simulation. See the simulation result in the below interactive Fig. A. 5 
to understand the saddle geometry that results from bending the Miura-ori. Readers may 
want to play with different toolbar options to better visualize the geometry. 
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FIG. A. 5: 3D geometry of a bent Miura plate made of 21 by 21 unit cells with a = 6 = 7r/3. 
For better display purpose, we use an example with pronounced deformation. However in the 
simulation we have done, we make sure that the radius of curvature is at least 10 times larger than 
the plate size, such that the deformation is within linear regime. Readers may want to play with 
different toolbar options to better visualize the geometry. 
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